Pinning/depinning of crack fronts in heterogeneous materials 
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The fatigue fracture surfaces of a metallic alloy, and the stress corrosion fracture surfaces of glass 
are investigated as a function of crack velocity. It is shown that in both cases, there are two fracture 
regimes, which have a well defined self-affine signature. At high enough length scales, the universal 
roughness index £ ~ 0.78 is recovered. At smaller length scales, the roughness exponent is close to 
£ c ~ 0.50. The crossover length £ c separating these two regimes strongly depends on the material, 
and exhibits a power-law decrease with the measured crack velocity £ c oc v~^ ', with tj> ~ 1. The 
exponents v and f3 characterising the dependence of £ c and v upon the pulling force are shown to 
be close to v ~ 2 and /3 — 2. 
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The pinning/depinning transition [gjj2j has been the 
subject of many theoretical studies in the recent years 
H . But although there are in principle numerous ap- 
plications, experimental examples remain scarce. In this 
letter, we analyse experimental results concerning crack 
propagation in two very different materials — a metal- 
lic alloy and glass — within the framework of models of 
pinning/depinning of moving lines through randomly dis- 
tributed obstacles • 

Fracture surfaces of many heterogeneous materials 
have been studied, with different experimental tech- 
niques. As shown first by Mandelbrot et al M, these 
surfaces are self-affine Q, with a roughness index £ in 
most cases close to the value 0.8. This quantity was later 
conjectured to be universal |^,^|, i.e. independant of the 
material and of the fracture mode. As far as metallic 
materials are concerned, this exponent has been recently 
characterized over five decades of length scales (0.5 nm- 
0.5 mm ) @PL This universality was first questioned by 
Milman ct al |12| on the basis of Scanning Tunneling Mi- 
croscopy experiments where fracture surfaces of metallic 
materials were investigated at the nanometer scale. The 
reported values of roughness exponents in the latter case 
were significantly smaller than 0.8, closer to 0.5. This has 
been interpreted jl0],[l3| as a kinetic effect similar to the 
one expected for a moving line near its depinning tran- 
sition. Indeed, it was proposed recently that a fracture 
surface could be modelled as the trace left by the crack 
front moving through randomly distributed microstruc- 
tural obstacles jlj]. 

Various models have been developped to calculate the 
critical exponents characterising the morphology of these 
lines either in 2 dimensions (moving front) or in 3 dimen- 
sions (moving front leaving a surface behind it). Erta§ 
and Kardar studied a local nonlinear three-dimensional 
Langevin equation to describe the morphology of poly- 
mers in shear flows or the motion of flux lines in super- 
conductors B, and it was conjectured that these models 



might also be relevant for fracture |14| . The line is pulled 
away with a constant force F. Non linearities account 
for the variations of the local crack speed with the local 
orientation of the front. This equation leads to a large 
number of regimes, depending on the relative values of 
the prefactors of the non linear terms. For some values, 
this model predicts that for a finite velocity v, the rough- 
ness exponent is 0.75 at "large length scales" and 0.5 at 
"short length scales", the two regimes being separated 
by a crossover length £ c . The short length scales regime 
corresponds to the vicinity of the depinning transition 
[GIB where the crack front is just able to free itself from 
the pinning microstructural obstacles. In this case, i.e. 
when F is close although higher than a critical force F c 
under which the line remains still, the velocity v tends 
to zero, v oc (F — F c )@ , and £ c diverges as £ c oc v^^. 
In this particular model 0], <f> = 3. It will be shown in 
the following that the observed behavior is in qualitative 
agreement with this scenario 10|, although the measured 
value of 4> is significantly smaller. 

Note that this behavior also corresponds to the re- 
sults of recent large-scale molecular dynamics simulations 
Jl5|-p7[ for amorphous materials. 

In this letter, quantitative results are presented, which 
show that £ c indeed decreases with the crack velocity, 
and lead to an estimate of <j>. Variations of the crossover 
length £ c with the average crack velocity are presented 
here both for the fatigue fracture of the TisAl-based 
Supera2 intcrmctallic alloy and for the stress corrosion 
fracture of soda-lime silica glass. 

Two notched compact tension specimens of Supera^ 
are broken in fatigue. Fatigue tests are carried out us- 
ing an electro-servohydraulic testing machine, operating 
under load control. The test is performed in air with a 
constant stress ratio R = a- min /a max — 0.1 (<r max and 
Urain are respectively the minimum and the maximum 
stresses), at a frequency / = 30Hz. The evolution of the 
crack length a with time is measured with the potential 
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drop method |18[ . The fracture surfaces are observed for 
4 different velocities spanning a wide range, using both 
an atomic force microscope (AFM) and a standard scan- 
ning electron microscope (SEM). The SEM observations 
consist in cutting and polishing the NiPd-plated frac- 
tured samples perpendicularly to the direction of crack 
propagation and registering images of the profiles at var- 
ious magnifications for each of the 4 regions (see [pjjf). 
For AFM observations, five profiles of length 1 /im are 
registered in each region. 
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Fig 1: Z max (r)j \f^ c is plotted against r/£ c for the 
two materials separately. Note that in these re- 
duced units, the plots corresponding to the various 
velocities collapse on the same curve. Although the 
crossover regions are quite different for the two ma- 
terials, the asymptotic regimes are well described by 
power laws with exponents 0.5 (r/£ c -C 1) and 0.78 
(r/& » !)• 

Fracture surfaces of soda-lime silica glass have been 
prepared by controlling the crack propagation with a 
four points bending system. After the initial propaga- 
tion, which allows to relax all residual stresses, the plate 
is properly loaded in order to obtain the required aver- 
age crack velocity. This velocity is measured by imaging 
the crack tip with AFM at different times. The humidity 
rate has been measured, and kept between 37 and 41%. 
The controlled crack propagation is maintained over a 
distance of about 30^m, so that fracture surfaces can 
be easily probed with AFM. Crack velocities range from 
2 10~ 9 to 10" 7 ms" 1 . Ten AFM height profiles of length 
1.5/im are registered perpendicularly to the direction of 
crack propagation, on three samples, and along this di- 
rection for four other specimens. As it will be shown in 
the following, no significant anisotropy could be detected. 

In order to determine the roughness exponents £ and 
the crossover length £ c of the profiles recorded, the 
Hurst method is used (2l|, where the following quan- 
tity is computed: Z max {r) = (max[z(r')] ro<r i <To+r — 
min[z(r')] ro<r i <ro + r } ro , Z max (r) oc r 1 *. r is the width of 



the window, and Z max (r) is the difference between the 
maximum and the minimum heights z within the win- 
dow, averaged over all possible origins r a of the window 
belonging to the profile. 




crack velocity (m/s) 



Fig 2: Evolution of the cross over length £ c with the 
crack velocity for the Supera^ (°) and for soda-lime 
silica glass (A/D: perpendicular/parallel to the di- 
rection of crack propagation). £ c is plotted versus v 
on a log-log plot, exhibitting a power-law dependence 
with an exponent <f> ~ 0.91. 

As far as the Supera2 is concerned, the simultaneous 
use of SEM and AFM allows for an observation of the 
fracture surfaces over 5 or 6 decades of length scales. 
Within the whole range of observations, Z max (r) is very 
well fitted by the sum of two power laws, Z max (r) = 
A(( r /£c)°' 5 + {r/£c) - 78 )- The small and large length 
scales roughness indices — 0.5 and 0.78 respectively — , 
are chosen to fit with the results of previous experiments 
JlOj l . On the contrary, the crossover between the two 
regimes is much sharper in the case of glass (which seems 
to be the case for other amorphous materials |19,[l7j), 
and in this case, the crossover length £ c is determined as 
the intersection of the two asymptotic power law regimes 
with exponents 0.5 and 0.78. Once the crossover lengths 
have been determined in each case, it is possible to plot 
Zmax as a function of r/£ c - In Fig. 1, the curves 
Z m ax{i") / relative to each material are plotted as a 
function of r/£ c and shown to collapse on the same mas- 
ter curve. In both cases, the asymptotic regimes are 
well described by power laws with exponents 0.5 at small 
length scales (r/£ c -C 1), and 0.78 at large length scales 
(r/£ c 3> 1). In other words, one can write: 

Z max {r) ex r - 5 /(f) (1) 

with f(x — > 0) ~ 1 and f(x » 1) ~ a; 28 , showing that 
the amplitude of the small length scales contribution is 
independent of crack velocity. 

The results obtained on materials as different as an 
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intermetallic alloy and a glass thus confirm previous ob- 
servations (l^Jl^l , where the short and large length scales 
regimes were interpreted, respectively, as a "quasi-static" 
and a "dynamic" regime. 

As it can be seen in Fig. 2, £ c decreases with the crack 
velocity v in both cases, although the measured values 
of £ c are approximately 1000 times larger in the case of 
the Supera^ than in the case of glass. The experimental 
results are compatible with a power law decrease in 
both cases. However, the estimated value of c/> is close to 
unity instead of 3. Note that the values of £ c measured 
for glass for the higher velocities might be overestimated: 
in this case the precision is very bad, since £ c is of the or- 
der of some nanometers, i.e. close to the limit resolution 
of the AFM. 



ity is not uniform, and intermittency is likely to occur. 
Thus, in both cases, the pinning/depinning scenario is 
qualitatively satisfactory. 

Fig. 4 shows the evolution of the crack velocity v as a 
function of AK — AKth for the Supera^ and as a func- 
tion of K — Kxh for glass. In the case of the Supera^, 
experimental measurements reveal a power law increase 
without any change between the so-called "threshold" 
and Paris regimes. On the contrary when static fracture 
occurs, a clear deviation from the power law can be ob- 
served for high values of AK — AKxh ■ A fit of these data 
gives j3 ~ 2. This value is compatible with the measure- 
ments on glass (Fig. 4). 
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Fig 3: Supertax : the fatigue crack velocity (□) as 
well as the cross over length £ c (o) are plotted against 
the stress intensity factor AK. The threshold value 
AKxh is indicated. 

In Fig. 3, both £ c and v are plotted against the stress 
intensity factor AK = (a m ax — Gmin)\/o- The fatigue 
crack growth in this regime is widely known as intermit- 
tent [p4U23] . In this regime, the crack tip opens and closes 
many times before it can extend over a small distance. 
This process is repeated several times and causes incre- 
mental crack advance. The number of cycles required 
to get the crack to advance decreases as AK increases, 
and the crack motion is more and more continuous, mi- 
crostructural obstacles being efficient at smaller length 
scales. At a given time, the force F exerted on the frac- 
ture front is proportional to AK, while the threshold 
force F c is proportional to AKxh (defined in Fig. 3). 
The frequency of oscillation of these forces being far more 
rapid than crack propagation, the average force only can 
be considered, which legitimates the analogy with the 
above-quoted models. In the case of glass, F is pro- 
portional to the stress intensity factor K, while F c is 
proportional to the threshold Kth- Preliminary results 
indicate that, in the sub-critical regime, the crack veloc- 
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Fig 4: Supera 2 (°): the fatigue crack velocity (white 
symbols) is plotted versus AK — AKtk on a log-log 
plot, as well as the crossover length £ c (black sym- 
bols). 

Glass (A/D perpendicular/parallel to the direction 
of crack propagation): the crack velocity is plotted 
as a function of K — Ktk (white symbols) , as well as 
£ c (black symbols). 

In Fig. 4, £ c is plotted also for both materials. A 
power law decrease can be observed, and the fit of the 
data relative to the metal gives v ~ 2, compatible with 
the results on glass. 

One can note that <f>= u//3 ~ 1. On the other hand, it 
is expected that the exponent n characterising the range 
of interactions (n = 2 for short range forces [JO and n = 1 
for long range ones f|) is related to the exponents v and 
£|| through the relation: 



c„ + i 



(2) 



where 0i is the in-plane roughness exponent |2q| . £\\ was 
determined previously on the Supera2 p6[ , and shown 
to be close to 0\ — 0.54. This leads to a value of 
n ~ 1.03 veryclose to unity as expected for elastic in- 
teractions fiE2,E3| . Note that the 2d model of Thomas 
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and Paczuski |^3| leads to £ = 0.5 and v = 2, but also to 

Knowing /3 and v, one can in principle deduce the value 
of the dynamic exponents z» and z± describing the short 
time evolution of the front, repectively in the direction of 
crack propagation and perpendicularly to it. zh = C|| + § 
leads to zn ~ 1.5, while z± = zn + — should indicate 
that z± ~ 2. Hence, perturbations on the crack front are 
diffusive perpendicularly to the direction of crack propa- 
gation, while they are slightly hyper-diffusive along this 
direction. 

For ductile materials as the Supero^, the plastic zone 
size Rpiast should be a relevant length scale as well. Al- 
though the same behaviour is observed within the whole 
range of AKs, it can be noted that R p i as t overpasses £ c 
for the two experiments corresponding to higher veloci- 
ties. 

The transition between the threshold regime and the 
Paris regime might be associated respectively with £ c > 
Rpiast and £ c < R p i as t- It has been shown in Fig. 2 
that a coarser microstructure gives rise to a higher £ c for 
the same crack velocity. On the other hand, Yoder et al 
showed that, in the case of a titanium-based alloy, a 
coarser microstructure corresponds to a higher value of 
AK at which the transition occurs. This result supports 
the idea that the transition between the threshold regime 
and the Paris regime could be associated to a competi- 
tion between R p i as t and £ c . Note also that in the case 
of glass, the plastic zone size has been estimated [^7j to 
be of the order of some nanometers, i.e. of the order of 
magnitude of £ c . 

Further experiments on different materials are needed 
in order to confirm the general character of the pin- 
ning/depinning scenario, and to allow for a more pre- 
cise determination of the critical exponents. In order to 
investigate the role of plasticity, experiments will be per- 
formed on an aluminium alloy, for which plastic zones 
are much larger than in the case of the Supera2- 
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